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Low-power radio galaxy environments in the Subaru/XMM-Newton 
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ABSTRACT 

We present multi-object spectroscopy of galaxies in the immediate (Mpc-scale) environments 
of four low-power (ii.4GHz ^ 10^*^ WHz^^) radio galaxies at 2: ^ 0.5, selected from the 
Subaru/XMM-Newton Deep Field. We use the spectra to calculate velocity dispersions and 
central redshifts of the groups the radio galaxies inhabit, and combined with XMM-Newton 
(0.3-lOkeV) X-ray observations investigate the Lx-(Jv and Tx-<Jv scaling relationships. All 
the radio galaxies reside in moderately rich groups - intermediate environments between poor 
groups and rich clusters, with remarkably similar X-ray properties. We concentrate our dis- 
cussion on our best statistical example that we interpret as a low-power (FR I) source triggered 
within a sub-group, which in turn is interacting with a nearby group of galaxies, containing 
the bulk of the X-ray emission for the system - a basic scenario which can be compared to 
more powerful radio sources at both high (z > 4) and low (z < 0.1) redshifts. This suggests 
that galaxy-galaxy interactions triggered by group mergers may play an important role in the 
life-cycle of radio galaxies at all epochs and luminosities. 

Key words: cosmology; observations - galaxies: evolution - galaxies; radio 



' 1 INTRODUCTION 



In recent years, feedback mechanisms from radio galaxies (that is, 
galaxies which host radio-loud active galactic nuclei) have become 
the method of choice for curtailing the bright end of the galaxy 
luminosity function in models of galaxy formation, since models 
without feedback tend to over-produce the number of very lumi- 
nous galaxies compared to what is observed (e.g. Bower et al. 
2006). Many high-redshift radio galaxies lie in (proto-)cluster envi- 
rorraients (e.g. Venemans et al. 2003) and the energies provided by 
their radio jets (the bulk kinetic power of the radio jets is several or- 
ders of magnitude larger than the radio luminosity; e.g., Rawlings 
& Saunders 1991) is sufficient to shut down star forming activity 
throughout a cluster (Rawlings & Jarvis 2004). Furthermore, it has 
been suggested that powerful radio sources may be triggered by 
galaxy-galaxy interactions during the merging of subcluster units 
(Ellingson, Yee & Green 1991; Simpson & Rawlings 2002). 

Lower luminosity sources are also believed to be vital in shap- 
ing the galaxy luminosity function by providing a low power, but 
high duty cycle, mode of feedback (e.g., Croton et al. 2005; Bower 
et al. 2006; Best et al. 2006). This mode of feedback appears to be 
most common in the most massive galaxies, which are, of course, 
preferentially located in the centres of galaxy clusters. Therefore 
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there appears to be significant interplay between the formation and 
evolution of both the radio sources and their cluster environments. 
Studying the environments of low-power radio sources is likely 
a promising route towards understanding radio galaxy feedback. 
However, while the radio galaxy itself can influence the thermo- 
dynamic history of its environment (i.e. by depositing energy), the 
picture is complicated by global environmental processes such as 
the merging of sub-cluster units, which can temporarily boost the 
X-ray luminosity and temperature of the intracluster medium (e.g. 
Randall, Sarazin & Ricker 2002), making such systems more read- 
ily detectable in X-ray surveys (a matter of caution in using such 
X-ray detected systems as Cosmological probes). It is essential to 
perform detailed studies of the radio, optical, and X-ray properties 
of putative dense regions in the cosmic web to ascertain how merg- 
ers and radio source activity affect the life-cycle of clusters. 

The luminosity and temperature of a cluster are known to 
follow a scaling relationship (the Lx-Tx relation), Lx oc Tx, 
where a has been measured in the range ~2.7-3 (Edge & Stew- 
art 1991; David, Jones & Forman 1995; Allen & Fabian 1998; 
Markevitch 1998; Arnaud & Evrard 1999) This is at odds with 
that expected for clusters formed by gravitational structure forma- 
tion, with Lx oc Tx (Kaiser 1986). It is thought that the dispar- 
ity is due to the interruption of cooling by feedback from active 
galaxies, which can impart energy via outflow shocks (for example 
from radio-jets or super-winds), or by performing work by sub- 
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Table 1. Summary of radio galaxy targets. I.D.'s, positions and fluxes are taken from an early version of the 1.4 GHz Very Large Array radio mosaic presented 
in Simpson et al. 2006 (S06). In this paper for convenience we shortem lAU identifications of the radio galaxies to JEG 1^. We quote radio luminosities here 
for completeness, however we discuss the redshift determination in §3.2. 



I.D. 


lAU I.D. 




S06 1.D. 


R.A. 
(J2000.0) 


Dec. 
(J2000.0) 


cSXDF 
''l.4GHz 
(raJy) 


^1.4 GHz 
(1024 WHz-1) 


JEGl 


VLAJ021945 


-04535 




02 19 45.3 


-04 53 33 


11.86±0.13 


4.36±0.05 


JEG 2 


VLAJ021823 


-05250 


VLAOOll 


02 18 23.5 


-05 25 00 


7.95±0.10 


14.35±0.18 


JEG 3 


VLAJ021737 


-05134 


VLA0033 


02 17 37.2 


-05 13 28 


2.37±0.06 


4.26±0.11 


IEG4 


VLAJ021842 


-05328 


VLA0065 


02 18 42.1 


-05 32 51 


0.96±0.08 


0.48±0.04 
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Figure 1. /? — z' vs. z' Field-corrected colour-magnitude diagrams for each mask. Photometry was extracted from the SXDF catalog by choosing all galaxies 
within an aperture of radius 2.5' centered on the brightest target galaxy. Colours are measured in 2" apertures, and we also indicate a secondary colour, i' — z' 
(arbitrarily chosen), to emphasise the location of the red-sequence as larger points. The colours of targeted galaxies in this work are indicated by larger open 
circles. 



sonically inflating bubbles in the intracluster medium (ICM; Si- 
jacki & Springe! 2006). Best et al. (2005) suggest that, given the 
short life-times of radio-loud AGN (10-100 Myr), episodic trigger- 
ing must occur in order for feedback to quench cooling. However, 
much lower luminosity sources (Li.4ghz < 10^^ WHz~^) can 
impart energy into the IGM with a higher duty cycle than the more 
powerful, and rarer, sources (Best et al. 2006). 

McLure et al. (2004) have recently constructed a sample of 
radio galaxies at z ~ 0.5 which spans three orders of magni- 
tude in radio luminosity, but even the faintest of these sources 
are close to the Fanaroff-Riley break (Fanaroff & Riley 1974) 
and it is desirable to push to even fainter radio luminosities to 
fully investigate the radio galaxy-cluster symbiosis. In this pa- 
per we present multi-object spectroscopy of galaxies around four 
low-power (Z/i.4GHz < 10^^ WHz"'^) radio galaxies aX z ^ 0.5 
in the Subaru/XMM-Newton Deep Field (SXDF) (Sekiguchi et al. 
2001, 2006 in prep). The resultant spectroscopic redshifts allow 



us to investigate the environments of these radio sources which 
could provide a large proportion of the local volume averaged heat- 
ing rate if the AGN luminosity function is flat below Li.4ghz ~ 
lO^'' W Hz" ^ (Best et al. 2006). We combine these data with X-ray 
observations of the radio galaxies' environments in order to inves- 
tigate the global environmental properties - namely the Lx — Ov 
and Tx — o-^ relations. This parameter space can be used as a di- 
agnostic of the thermodynamic history of groups or clusters, since 
if an AGN is significantly feeding energy back into the intraclus- 
ter medium (ICM), then the group or cluster might be expected to 
deviate from the standard empirical Lx (Tx) — o-„ relation for mas- 
sive systems. Throughout we adopt a flat geometry with f2m = 0.3, 
= 0.7 and Ho = lOO/i kms"^ Mpc"^ where h = 0.75. All 
magnitudes in this paper are on the AB scale. 
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Figure 2. Successfully cross-correlated spectra for the radio galaxies field 
of (top-down) JEG 1^. The spectra are plotted in the rest-frame. As a guide, 
we identify as dotted vertical lines the spectral features [On], K, H, HS, G, 
H7, and H/3. We do not plot wavelengths blueward of 3500A or redward of 

5OOOA (restframe) since in these regimes the response and signal-to-noise 
of the spectra do not allow accurate analysis. Note that the shaded region 
higMghts a sky absorption feature. 



Table 2. Overview of the XMM-Newton observations. Tabulated are the 
Field IDs, the centres of the pointings and the exposiffe times of the rel- 
evant high-background-cleaned datasets for each camera: MOS 1, M0S2 
and PN. Note that the coordinates are in J 2000. 



Field 


R.A. 


Dec. 


Exp. time [MOS 1 MOS 2 PN] 




hh mm ss 


' " 


(ksec) 


SDS-la 


02 18 00 


-05 00 00 


40 


SDS-lb 


02 18 00 


-05 00 00 


42 


SDS-2 


02 19 36 


-05 00 00 


47 48 40 


SDS-6 


02 17 12 


-05 20 47 


49 49 47 


SDS-7 


02 18 48 


-05 20 47 


40 41 35 



sequence were selected for spectroscopic follow-up. This selection 
was simply performed 'by eye' - if the colour-magnitude diagram 
displayed an overdensity of points at 18 < mna < 20.5 with tight 
scatter in colour space. Note that local X-ray properties were not 
taken into account when choosing the target radio galaxies - this is 
simply an optical selection of radio galaxies. 

Masks were designed using the LDSS-2 mask preparation 
software, with objects along the red sequence (typically using a 
colour cut of width 0.3 magnitudes) included in the input file (Fig- 
ure 1). This software allows relative priorities to be assigned to tar- 
gets, which were set to be equal to the magnitudes of the galaxies, 
to increase the likelihood of successfully determining redshift. 

Observations were conducted using the medium-red (300 
lines/mm, Ao = 5550A) grism of LDSS-2 in multi-object mode 
on the 6.5-m Magellan telescope in Las Campanas, Chile over the 
nights 17-20 October 2003 (effective resolution 13.3A). The MOS 
frames were trimmed for overscan, flat-fielded, rectified, cleaned 
and wavelength calibrated using a sequence of Python routines (D. 
Kelson, priv. comm.); sky-subtraction and spectrum extraction was 
performed in IDL. Due to a lack of standard-star observations, for 
flux-calibration and to correct for the wavelength dependent re- 
sponse of the detector, we extracted BVRi'z photometry for the 
targets from the Subaru/XMM-Newton Deep Field (SXDF) catalog. 
Using the five sources in each mask with the strongest signal-to- 
noise ratio (SNR), we fit mask-specific response curves to the 1-D 
spectra, which are then applied to each extracted target spectrum. 
Although flux-calibration and response correction is not necessary 
for redshift determination via cross-correlation, we perform these 
steps in order to present a sample of the spectra in Figure 2. 



2 OBSERVATIONS & REDUCTION 
2.1 Sample selection and spectroscopy 

Targets were selected from an early version of the 1.4-GHz Very 
Large Array radio mosaic presented in Simpson et al. (2006), and 
listed in Table 1. The limited red sensitivity of the Low Disper- 
sion Survey Spectrograph (LDSS-2, AUington-Smith et al. 1994) 
placed an upper limit of 2: ^ 0.7 as appropriate for our study, so 
we selected radio sources whose optical counterparts are resolved 
and have 19 < i? < 20.5, since the tightness of the Hubble di- 
agram for radio galaxies allows redshifts to be estimated from a 
single magnitude, and we aim to target moderate redshift environ- 
ments (Cruz et al. 2007). Colour-magnitude diagrams were pro- 
duced for 2.5' regions around each radio galaxy and the four opti- 
cally brightest radio galaxies whose environments displayed a red 



2.2 X-ray data 

The SXDF incorporates a deep, large-area X-ray mosaic with 
XMM-Newton, consisting of 7 overlapping pointings covering 
^l.Sdeg'^ region of the high-Galactic latitude sky with an expo- 
sure time of looks in the central field (in separate exposures) and 
50 ks in the flanking fields (see Table 2). Four of the pointings were 
carried out in August 2000, and the remaining three were made in 
August 2002 and January 2003. 

All the XMM-Newton EPIC data, i.e. the data from the two 
MOS cameras and the single PN camera, were taken in full-frame 
mode with the thin filter in place. These data have been reprocessed 
using the standard procedures in the most recent XMM-Newton 
SAS (Science Analysis System) — v.6.5. All the datasets where 
at least one of the radio galaxies was seen to fall within the field- 
of-view (FOV) of the detectors were then analysed further. JEG 3 
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was observed in pointings SDS-la (PN), SDS-lb (PN) and SDS- 
6 (Ml, M2, PN) (part of the field lies off the edge of the MOS 
FOVs in SDS 1 a and SDS lb). Both JEG 2 and JEG 4 were observed 
in pointing SDS-7 (MOS 1, MOS 2, PN). Finally JEGl was ob- 
served in pointing SDS-2 (MOS 1, MOS 2, PN). Periods of high- 
background were filtered out of each dataset by creating a high- 
energy 10—15 keV Ughtcurve of single events over the entire field 
of view, and selecting times when this Ughtcurve peaked above 
0.75 ct s"^ (for PN) or 0.25 ct s"^ (for MOS). It is known that soft 
proton flares can affect the background at low energies at different 
levels and times than high energies (e.g. Nevalainen et al. 2005). 
This was not seen to be the case in the datasets analysed here how- 
ever, as full energy-band (1-10 keV) lightcurves (extracted from 
off-source regions) showed essentially identical times of flaring ac- 
tivity as the high energy lightcurves. The high-background times in 
each dataset are seen to be very distinct and essentially unchang- 
ing with energy-band. The exposure times of the resulting cleaned 
relevant datasets are given in Tabic 2. 

Source spectra were extracted from the relevant datasets from 
apertures centred on the source positions. Extraction radii, esti- 
mated from where the radial surface brightness profiles were seen 
to fall to the general surrounding background level, were set to 50", 
21", 36" and 47" for JEG 1^ respectively. Note that the surface 
brightness profiles were co-added for all 3 EPIC cameras for im- 
proved S/N. Standard SAS source-detection algorithms, uncovered 
a large number of X-ray sources, both point-like and extended. Data 
due to contaminating sources lying within the extraction regions 
needed to be removed from the spectra. This proved necessary only 
for JEG 3, where two sources were detected just within the extrac- 
tion region. Data from the two sources were removed to a radius 
of 28" (and the extraction area is thus reduced by ~ 10%). Back- 
ground spectra were extracted from each cleaned dataset from an 
annulus of inner radius 60" and outer radius 180" around each po- 
sition. As point sources were seen to contaminate these larger-area 
background spectra, the data from these detected sources were re- 
moved from the background spectra to a radius of 60". We have 
verified that all the X-ray sources are clearly extended using both 
simple (sliding box) and more sophisticated wavelet analysis tech- 
niques. The possibility that JEG 2 is an unresolved source com- 
plex cannot of course be completely ruled out given its low signal 
to noise. Given the relative weakness of all of the X-ray sources, 
we consider that the extraction radii used are close to optimvun to 
provide the best signal to noise in the extracted spectra. Ancillary 
Response Files (ARFs) containing the telescope effective area, fil- 
ter transmission and quantum efficiency curves were created for the 
cluster spectra (taking into account the fact that the extraction re- 
gions are extended), and were checked to confirm that the correct 
extraction area calculations (complicated with the exclusion of con- 
taminating point sources, etc.) had been performed. Finally Redis- 
tribution Matrix Files (RMFs) files were generated. These files de- 
scribe, for a given incident X-ray photon energy, the observed pho- 
ton energy distribution over the instrvunent channels. These can be 
used in conjunction with the ARFs to perform X-ray spectral anal- 
ysis - we use the most commonly used analysis package: XSPEC. 
The source spectral channels were binned together to give a mini- 
mum of 20 net counts per bin. A detailed discussion of the X-ray 
spectral fitting and the results is given in §3.4 
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Figure 3. Relative velocity distribution of galaxies where successful 
redshift-determination was performed. Relative velocities are in the cluster 
frame, over a range ±2000 km s~ ^ . Insets show the distribution in redshift 

space, over a wider velocity range. We overlay Gaussian profiles based on 
the fitted velocity dispersion in each case. 



3 ANALYSIS 



3.1 Colour-magnitude diagrams 

In Figure 1 we plot the R — z' vs. z' colour-magnitude rela- 
tions (CMRs) for each field, which have been constructed from 
the photometry of all galaxies within 2.5' of each radio-galaxy tar- 
get. We correct the colour-magnitude diagrams for field contam- 
ination using a method similar to Kodama et al. (2004) and Ko- 
dama & Bower (2001), which we will briefly describe here. First, 
colour-magnitude space is divided into a coarse grid (Az' = 2, 
A{R — z') = 0.8), and the number of galaxies within each grid- 
square counted: A'^^''' . An identical count is performed for a control 
field, Aff/"'"°', where control galaxies are selected within an aper- 
ture of identical radius to the one used for each target mask, but cen- 
tered on a random point on the sky (but within the catalog). We then 



ics within it are field members: f 



^controY(^gal_ 



rcontrol \ 

For each galaxy in colour-magnitude space, we pick a random num- 
ber, p, in the interval {0,1} and compare it to Pij-, a galaxy is re- 
moved from the diagram if p < Pij. In order to reduce effects 
of cosmic-variance (e.g. randomly picking another cluster or void 
as the control sample), we build up an average probability map by 
repeating the control sample count and probability calculation 10 
times around different random positions in the catalog. 

There are clear red-sequences around all four radio galaxies, 
and to emphasise this, we also plot a secondary colour i — z', and 
the locations of galaxies in colour space targeted for spectroscopy. 
The secondary colour-cut is chosen merely to highlight the red- 
sequence. 
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Figure 4. The XMM-Newton EPIC spectra of (row-wise from top left) JEG 1^, together with the best-fitting thermal plasma model fitted using XSPEC. 
differences are shown in the lower panels of each plot. The colour-codings correspond to the MOS 1, MOS 2 and PN cameras and field pointings (see Table 2 
and §2.2). For JEG 1 (SDS-2), JEG2 (SDS-7) & JEG4 (SDS-7) black is MOS 1, red is MOS 2 and green is PN. For JEG 3 black is PN (SDS-la), red is PN 
(SDS-lb), green is MOS 1 (SDS-6), blue is MOS 2 (SDS-6) and cyan is PS (SDS-6). The resulting cluster parameters of Lx and Tx are presented in Table 4. 



3.2 Redshift determination and velocity dispersions 

For redshift determination we use the VELOCITY code of Kelson 
et al. (D. Kelson, priv. comm.), which cross-correlates target and 
template spectra. The errors are determined from the topology of 
the cross-correlation peak (i.e. the width of a Gaussian profile fit 
to the peak in the cross-correlation spectrum) and the noise associ- 
ated with the spectrum; typically in the range 30-210kms^^. We 
present spectra of the target radio galaxies which have been suc- 
cessfully cross-correlated and shifted to the rest-frame in Figure 2. 
In Table 3 we list the coordinates, z' magnitudes and redshifts of 
successfully cross-correlated targets, with the radio galaxies high- 
lighted. 

To calculate the velocity dispersions of each cluster, we fol- 
low the iterative method used by Lubin, Oke & Postman (2002) 
and others (e.g. Willis et al. 2005). Initially we estimate that the 
cluster central redshift is that of the target radio galaxies, and se- 
lect all other galaxies with \Az\ < 0.06 in redshift space. We 
then calculate the biweight mean and scale of the velocity distri- 
bution (Beers, Flynn & Gebhardt 1990) which correspond to the 
central velocity location, Vc, and dispersion, av of the cluster. We 
can use this to calculate the relative radial velocities in the cluster 
frame: Av = c{z — Zc)/(1 + Zc). The original distribution is re- 
vised, and any galaxy that lies >3(t„ away from Vc, or has | Aii| > 
3500 km s~^ is rejected from the sample and the statistics are re- 



calculated. The final result is achieved when no more rejection is 
necessary, and the velocity dispersion and cluster redshift can be 
extracted. The results are presented in Table 4, with 1-cr errors on 
the cluster redshift and dispersion corresponding to 1000 bootstrap 
resamples of the redshift distribution. 

In Figure 3 we plot the redshift distribution of those galaxies 
for which we performed successful cross-correlations. We also plot 
the velocities of those galaxies within \ Av\ < 2000 km s^^ of the 
cluster central redshift, overlaid with standard Gaussian profiles. 
JEG 3 is the largest sample of cluster galaxies (in part due to the 
large number of galaxies targeted in this field). JEG 1 and JEG 4 
(at similar redshifts) have relatively few members after the process 
of rejection described above. JEG 2 has a high velocity dispersion, 
indicating a non-relaxed system, despite the fairly prominent red- 
sequence in this region. In part, the high velocity dispersion might 
be attributed to an outlying galaxy in the velocity distribution, but 
as we discuss in §3.3, JEG 2 might be part of a larger, more com- 
plicated stucture than a simple group. 

3.3 Environmental richness estimates 

Although a rigorous spectroscopic environmental analysis of the 
clusters is not feasible with the sample sizes in this study, we can 
make an estimate of the richness of the clusters by some sim- 
ple counting statistics using the photometric information from the 
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Table 3. Spectroscopic targets with total z' band magnitudes and cross-correlated redshifts. Redshift uncertainties are derived from the topology of the 
cross-correlation peak. The radio galaxy is highlighted in bold type in each case. 



I.D. 



RA 
(J2000.0) 



Dec. 
(J2000.0) 



Redshift 



Member* 



(mag) 













JEGl 












SXDF-iE-085449 


02 


19 


45 


252 


-04 


53 


33 


08 


18 


27 


SXDF-iE-079345 


02 


19 


39 


461 


-04 


53 


18 


56 


18 


83 


SXDF-iE-073826 


02 


19 


36 


864 


-04 


53 


29 


98 


18 


01 


SXDF-iE-094234 


02 


19 


48 


300 


-04 


53 


23 


65 


19 


56 


SXDF-iE-089522 


02 


19 


46 


118 


-04 


51 


22 


25 


19 


36 


SXDF-iE-085293 


02 


19 


43 


190 


-04 


52 


43 


50 


19 


33 


SXDF-iE-067505 


02 


19 


33 


689 


-04 


51 


45 


47 


18 


35 


SXDF-iE-080550 


02 


19 


41 


508 


-04 


52 


30 


60 


17 


61 












JEGl 










SXDF-iS-087629 


02 


18 


23 


532 


-05 


25 


00 


65 


19 


17 


SXDF-iS-105254 


02 


18 


19 


109 


-05 


23 


06 


32 


20 


02 


SXDF-iS-10()288 


02 


18 


19 


066 


-05 


23 


45 


74 


20 


19 


SXDF-iS-111951 


02 


18 


20 


074 


-05 


22 


14 


13 


21 


49 


SXDF-iS-108679 


02 


18 


20 


374 


-05 


22 


36 


40 


20 


64 


SXDF-iS-087515 


02 


18 


23 


777 


-05 


25 


49 


89 


21 


73 


SXDF-iS-()9()262 


02 


18 


26 


772 


-05 


25 


14 


35 


19 


66 


SXDF-iS-()7()840 


02 


18 


27 


202 


-05 


27 


36 


30 


19 


70 


SXDF-iS- 102854 


02 


18 


29 


954 


-05 


23 


30 


49 


20 


73 












JEG3 










SXDF-iS-170569 


02 


17 


37 


193 


-05 


13 


29 


61 


19 


30 


SXDF-iS-157505 


02 


17 


24 


230 


-05 


15 


19 


53 


19 


76 


SXDF-iS-173222 


02 


17 


26 


530 


-05 


13 


45 


93 


21 


01 


SXDF-iS-182658 


02 


17 


27 


502 


-05 


11 


45 


34 


21 


44 


SXDF-iS-168040 


02 


17 


29 


052 


-05 


12 


59 


61 


19 


89 


SXDF-iS-154410 


02 


17 


33 


286 


-05 


15 


51 


00 


19 


46 


SXDF-iS- 184999 


02 


17 


31 


322 


-05 


12 


17 


11 


20 


85 


SXDF-iS-172233 


02 


17 


35 


309 


-05 


13 


30 


71 


19 


76 


SXDF-iS-172341 


02 


17 


43 


358 


-05 


13 


30 


69 


20 


78 


SXDF-iS-172053 


02 


17 


34 


202 


-05 


13 


39 


44 


19 


36 


SXDF-iS-181772 


02 


17 


25 


450 


-05 


11 


54 


43 


19 


81 


SXDF-iS-169690 
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1.847, see Figure 7 and §4.2 
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Table 4. Spectroscopic, environmental and X-ray properties of target clusters. The cluster redshift Zc and line-of-sight velocity dispersion (T„ are calculated 
using the biweight mean and scale of the distribution (§3.2), with Icr errors corresponding to 1000 bootstrap resamples of the redshift distribution. The 
environmental richness statistics (see §3.3) - are the A^q.s counting statistic and Bgc, the amplitude of the galaxy-cluster cross-coiTelation function. We also 
tabulate the results of fitting thermal (mekal) plasma models to the X-ray spectra of the clusters (see §3.4 for details). The uncertainties represent 90% for a 
single interesting parameter, determined during the spectral fitting with the XSPEC software. 



Redshift properties Environmental properties X-ray properties (0.3-10 keV) 



I.D. 


2c 


(kms-1) 


fa. 

J mem 


A'o.s 


(Mpcl-'''^) 


kT 
(keV) 


Z 

(Zq) 


Lx 
(10^6 -^-j 




Net counts 


JEGl 
JEG2 
JEG3 
JEG4 


0.333±0.009 
0.649±0.001 
0.648±0.001 
0.38±0.01 


643±223 
1042±394 
774±170 
436±89 


5/8 
7/9 
13/18 
3/11 


24 ±5 
4±2 
23 ±5 
23 ±5 


599 ± 198 
696 ± 557 
1782 ± 585 
388 ± 204 


1 ^9+0.65 

'■^^-0.84 
, 07+1.00 
'■^'-0.59 

1 06+'^-52 
l.UO_Q 27 


^0.2 
uncon. 
^0.7 
;S0.8 


(-2+0.07 
"•"^-0.07 

, 7Q+0.29 
'•'^-0.83 
(,^,+0.06 
^■°-'-0.11 


57.6/64 

1.9/3 
51.2/51 
45.2/43 


477.29 ± 37.68 
43.16 ± 12.12 
328.21 ±33.51 
232.77 ±30.03 



^ fraction of members associated with the radio galaxy environment for successful redshifts 



SXDF catalog. The A'o.s measure is perhaps the best tool for a first- 
order quantitative analysis of a cluster environment - it only relies 
on an approximate cluster redshift (in order to calculate the angular 
radius of the counting aperture), and a satisfactorily deep photomet- 
ric catalog (Hill & Lilly 1991). It has been shown to correlate well 
with the more complex _Bgg statistic (Longair & Seldner 1979), 
which is a more rigorous measure of the density of galaxies in space 
when spectroscopic data is lacking, but when the shape of the lumi- 
nosity function is known (see uses in e.g. Farrah et al. 2004, Wold 
et al. 2000, 2001). In physical terms _Bgg is the amplitude of the 
galaxy-galaxy spatial correlation function: ^(r) — Bggr"^, where 
7 is typically 1.8-2.4 (Lilje & Efstathiou 1988, Moore et al. 1994, 
Croft, Dalton & Efstathiou 1999), though the low end of the range 
is typically chosen, with 7 — 1.77 (e.g. Yee & Ellingson 2003). 
Both these measures can be translated to the more familiar Abell 
richness scale. 

is calculated by choosing a target and counting the num- 
ber of galaxies within 0.5 Mpc in the magnitude range (m, m + 3), 
where m is the magnitude of the target galaxy, which in this case 
we choose to be the radio galaxies. A field correction is applied 
by subtracting the expected number of galaxies in this magnitude 
range from a control field (this control field was different for each 
target, chosen to be centred on a random point, > 1 Mpc from the 
radio galaxy). Note that the N0.5 statistic will be inaccurate in 
cases where there is strong differential evolution between the radio 
galaxy luminosity and companion cluster galaxies, but this can be 
compensated for by choosing the second or third brightest galaxy 
in the aperture and applying the same method. More importantly, it 
is clear from the spectroscopic results that the A'o.s statistic is con- 
taminated by line-of-sight structures not associated with the phys- 
ical environment of the radio galaxies. This is particularly evident 
for JEG 4 - the slightly higher redshift structure suggested by the 
number of z ~ 0.4 galaxies have magnitudes and colours that not 
only re-inforce the original red-sequence selection, but will be in- 
cluded in the A'^0.5 calculation. Thus it is important to note that this 
method can only provide a very rudimentary environmental analy- 
sis - a more sophisticated method is required. 

For this reason, we also calculate Bgc - where 'gc' stands for 
'galaxy-cluster centre' . This is a potentially useful statistic: derived 
solely from photometry, it has been shown to be in excellent agree- 
ment with its spectroscopic counterpart, and can be used as a pre- 
dictor for global cluster properties such as the velocity dispersion, 
virial mass and X-ray temperature (Yee & Ellingson 2003). Thus, 
we have another parameter to compare with our spectroscopically 
derived cluster velocity dispersions and X-ray observations, as well 



as a quantitative description of the richness of the radio galaxies' 
environments. 

A full derivation of the _Bgc statistic was made by Longair & 
Seldner (1979), and has been used in a practical sense by several 
other authors (see Wold et al. 2000, 2001, Farrah et al. 2003, Yee & 
Lopez-Cruz 1999), so we will not give a detailed procedure here, 
but the form of the statistic is given by: 



where pg is the surface density of field galaxies brighter than a limit 
T7^lim, de is the angular diameter distance to z, is an integration 
constant (7^=1.77 ~ 3.78) and ^(miini, z) is the integral luminos- 
ity function to a luminosity corresponding to mum at z. Age is the 
amplitude of the angular correlation function, calculated by com- 
paring the net excess of galaxies within 0.5 Mpc of the target with 
the number of background sources expected for an identical aper- 
ture in the field: Age = (iVnot/iVbg)((3 - 7)/2)6l^~\ The con- 
version to Bgc effectively deprojects from the celestial sphere 
into 3-D space. 

Perhaps the most important factor to consider is the choice of 
luminosity function, although Bgc is reasonably tolerant of incor- 
rect parameters (as long as they are within ~ 20% of their true val- 
ues). We estimate the shape of the luminosity function at z — 0.35 
and z = 0.65 by using the semi-analytic catalog output of the Mil- 
lenium SimulatiorQ (Bower et al. 2006, Springel et al. 2000). At 
each epoch, we fit a Schechter function to the absolute (simulated) 
_R-band magnitude counts, fixing the faint end slope a = —0.9. 
The resulting {AU, (p^,) are (-20.87 mag, 0.0042 Mpc"^) and 
(-20.65 mag, 0.0021 Mpc"^) at z = 0.35 and z = 0.65 respec- 
tively. The relevant luminosity function used in equation 1 is eval- 
uated from A/* + 3 to Mi, — 3, where the faint end is in a relatively 
flat part of the distribution, and the bright end is expected to include 
most of the cluster galaxies (note that changing these limits by ±1 
magnitude does not change the value of Bgc by more than la). The 
results are presented in Table 4. 

It is important to note that JEG 2 appears to have an environ- 
ment almost consistent with the field (although, as noted the No. 5 
statistic is likely to suffer contamination). This is inconsistent with 
the apparent red-sequence around the radio galaxy, which would 
suggest a moderate overdensity. However, we have seen that the 
very high velocity dispersion for this group, and structure in the ve- 
locity distribution might affect the richness estimates if the group 

|http : / / galaxy-catalogue . dur . ac . u"k| 
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is spatially spread out. JEG 1 and JEG 4 are moderately rich en- 
vironments (they have very similar levels of overdensity, although 
note that JEG 4 is contaminated by a slightly higher redshift line- 
of-sight structure, and so its true richness is likely to be slightly 
lower),consistent with rich groups. JEG 3 is also an overdensity, 
but the statistics suggest that the cluster is no more rich than JEG 1 
or JEG 4. Thus, the radio galaxy selection is detecting moderately 
rich groups of galaxies. 



3.4 X-ray spectral fitting 

The X-ray spectral fitting package xspec was used to fit X-ray 
spectral models to the cluster spectra. In terms of X-rays from 
galaxy clusters, it is commonly believed that the emission is pre- 
dominantly due to the hot cluster gas trapped in the potential well, 
absorbed by the hydrogen column density within our own Galaxy 
(i.e. little in the way of internal cluster absorption is expected). The 
spectra therefore were fitted with a standard model comprised of a 
mekal hot plasma emission spectrum, together with a wabs photo- 
electric absorption model (Morrison & McCammon 1983). Due to 
the small number of counts, the hydrogen column of the absorp- 
tion was fixed at the Galactic hydrogen colutim density in the line 
of sight to the Subaru pointings: Nu = 2.55 x 10^" atoms cm^^ 
(Dickey & Lockman 1990). Also the redshifts were fixed at the val- 
ues given in Table 4. The temperature and metallicity were left free 
to optimize. Acceptable spectral fits were found for all the clusters 
(though only a small number of counts is seen from JEG 2), and 
these results are given in Table 4. We present the best fit gas tem- 
perature and metallicity (with respect to the Solar value), the errors 
being 90% for a single interesting parameter, and the (0. 3 — 1 keV) 
absorption-corrected X-ray luminosity. Also tabulated are the best 
fit values and the number of degrees of freedom. For none of 
the clusters were there sufficient counts to warrant the use of a 
more complex model. We attempted to rebin the spectra into five 
counts per bin, with a fit based on Cash statistics (i.e. not subtract- 
ing the background in order to maintain counting statistics. Cash 
1979), however this does not provide stronger constraints on T or 
Z compared to the adopted method. The data and best-fit models 
are shown in Figure 4. 



4 DISCUSSION 

With the exception of JEG 2, the radio galaxies in each of our 
fields have environmental properties consistent with moderately 
rich groups at z ~ 0.5, in terms of their density statistic and X- 
ray properties. Despite the presence of a clear red sequence in the 
field of JEG 2 (Fig. 1), a simple A'o.s environmental measure sug- 
gests that the environment is consistent with the field. However, it 
is likely that this is in fact a cluster of similar richness to JEG 1, 
JEG 3 and JEG 4 given that there are several nearby galaxies with 
colours and spectroscopic redshifts close to the radio galaxy. The 
low-number statistics for this system makes it difficult to make any 
robust comments about the radio-galaxy environment, except that 
this is probably a dynamically young group. There is some evi- 
dence for a sheet at z ~ 0.65 (Simpson et al. in prep) and JEG 2 
may simply lie in an overdensity within this stucture. A large pro- 
jected distribution of red galaxies aroimd this group would in part 
explain the failure of the Bgo statistic to detect an overdensity: it 



could be possible that the control aperture itself is contaminated by 
the larger stucture. 

Our observations of JEGl and JEG 4 also do not contain 
enough confirmed spectroscopic members at the radio galaxies' 
redshifts to comment on their environments from a spectroscopic 
standpoint. The wide-angle radio tail of JEG 1 suggests rapid move- 
ment through a dense mediimi (806), consistent with a group-group 
merger (e.g. Jetha et al. 2006). There is insufficient spectroscopic 
data to confirm this though. Furthermore, in the case of JEG 4, our 
results could be affected by a background group, although it is not 
clear that this is contributing significantly to the X-ray Ivuninos- 
ity. This slightly more distant group could also be re-enforcing the 
strength of the red sequence (Fig. 1), and is a cautionary point re- 
garding selection in this way. 

4.1 Scaling relations 

The X-ray luminosities derived by fitting thermal models to the X- 
ray spectra are remarkably constant over our radio galaxy sample, 
with the environment of JEG 3 the hottest and intrinsically the most 
luminous (we note that the small number of counts seen from JEG 2 
results in a poorly constrained model). In Figure 5 we compare 
these clusters' properties with the Lx-Tx relation of Helsdon & 
Ponman's (2000) study of galaxy groups and low-luminosity clus- 
ters, more luminous and massive systems from Wu, Xue & Fang 
(1999), as well as more intermediate redshift (z < 0.6) X-ray se- 
lected groups/clusters from the XMM-Latge Scale Survey (XMM- 
LSS; WiUis et al. 2006). 

Firstly, the four clusters all have high-Lx and high-T com- 
pared to the group sample, indicating that they bridge the gap be- 
tween groups and rich clusters - a relatively unexplored region of 
parameter space. This is likely due to selection effects, thus there 
is potential in investigating these types of environment using a se- 
lection technique similar to the one outlined in this work. However, 
we note that the radio selected environments are quite similar to 
X-ray selected groups/clusters at similar redshifts, detected in the 
XMM-LSS (WiUis et al. 2006). We show also the disparity between 
the slopes of the Lx-Tx relation for groups (Helsdon & Ponman 
2000) and rich clusters (Wu, Xue & Fang 1999), which differ by 
a factor r^2. This is probably a real effect - the steepening of the 
slope for groups reflects the fact that these systems have a propsen- 
ity to exhibit deviations from the relation seen for the more massive 
systems due to processes other than cooling (Willis et al. 2005; 
Ponman, Cannon & Navarro 1999). Note that all the empirical re- 
lations found for groups and clusters are steeper than that expected 
for pure gravitational structure formation (Kaiser 1986); this im- 
plies that feedback as a source of energy input is important over all 
scales of clusters, but is more dominant in the low mass systems. 

One can also compare the cluster X-ray luminosity and ve- 
locity dispersion with the Lx-cr relations for these environments 
(Fig. 5). Although our small numer statistics provide only loose 
constraints on our environmental parameters, all our environments 
appear to have slightly large velocity dispersions compared to the 
temperature of the intracluster gas, although this is not statisti- 
cally significant. Once again, the velocity dispersions and temper- 
atures reveal that these radio galaxies inhabit intermediate environ- 
ments between groups and clusters. A similar picture is revealed 
by the Lx-ffv relationship, with all clusters having slightly large 
(T„ for their Lx- This is in good agreement with the observed ve- 
locity dispersions for other radio/optically selected samples (e.g. 
Rasimiussen et al. 2006; Popesso et al. 2007). 

For the remainder of the discussion we concentrate on JEG 3, 
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Figure 5. The correlation of global cluster properties: (top left) Lx-Tx (bottom left) Tx-^v (top right) Lx-Ov- We compare our results with both rich 
clusters (Wu, Xue & Fang 1999) and poor groups (both the 'loose' and 'compact' sample from Helsdon & Ponman 2000), as well as more recent results from 
the XMM-Large Scale Structure (ZMM-LSS) survey (Willis et al. 2005). In all parameter spaces shown, the sample in this work occupy the intermediate space 
between the rich/poor sample and are similar to the XMM-LSS sample. The dotted lines are the relationships derived for the groups sample from Helsdon & 
Ponman (2000), and the dashed lines show the equivalent relations derived using the rich cluster data from (Wu, Xue & Fang 1999). For clarity, we do not plot 
eiTors on the literature points (but these are taken into account in the linear fits). The X-ray data from this work are corrected for aperture losses. To estimate 
this, we extrapolate a King surface brightness profile of the form: S (r) = 5o [1 + {t Jtcoyg )2] -3/3+0.5 Assuming /3 = 0.5 and ^'corc — 200 kpc, we find that 
the apertures used recover ~80% of the total flux, extrapolating to 2r'corc. AH points are coirected to our adopted cosmology, with Hq = 75 km s~^ Mpc~^ . 



which has a spectroscopically-confirmed sample large enough to 
comment on the Mpc-scale environment of the radio galaxy. 



4.2 Discussion of JEG 3 

Optically, the radio galaxy is intriguing: it appears to be morpho- 
logically disturbed, with a blue and red part, suggestive of a recent 
or ongoing interaction triggering star-formation. In fact, closer in- 
spection of the spectra (and images) provide unambiguous evidence 
that the blue compenent is actually a strongly lensed background 
AGN or starburst galaxy al z = 1.847 - not associated with the 
group (Figure 7). This high redshift object does not affect our re- 



sults in any way, other than the fact that without the spectroscopic 
confirmation, this might be classified as a merger from the opti- 
cal images alone. In terms of the radio emission, we are confident 
that the majority of the emission originates in the lensing source - 
i.e. the original selection of this target is still secure (the morphol- 
ogy and position of the radio emission support this, see Fig. 6). In 
terms of the impact from X-ray emission from the lensed galaxy, 
we also estimate that our measurements of the X-ray luminosity of 
the group will not be affected. The main bulk of the X-ray emis- 
sion appears to be offset from the radio galaxy itself (Fig. 6) and 
the data were cleaned for point source emission prior to spectral 
fitting (§2.2). In the case that the point source is not cleaned, we es- 
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Figure 6. Here we present true colour image constructed from Subaru BVi-hand 8' X 8' images centered on the target radio galaxies JEG 1^ (orientated 
North up and East left). In addition to the colour information, we overplot as contours the radio galaxy morphology from the VLA images (white solid), and 
X-ray (0.3-10 keV) contours (green dashed). The radio contour levels have been chosen to emphasise the lobe morphology (note that these are VLA B-aiTay 
maps with a beam size of 5' X 4" with P. A. = 170°, see S06), and the X-ray contours are slightly smoothed with a Gaussian kernel of width and at 
~2,3,5 & la above the background in the vicinity of the environment. Finally, using the redshift data, we highlight galaxies which have | Ad| < 2000 km s~^ 
relative to the cluster redshift. We plot positive values of At) as open circles and negative values as open squares. 



timate the impact on our measurement of Lx '■ if the lensed galaxy 
is a starburst, then we expect the X-ray luminosity to be of order 
~ 2 X 10^" W (based on an X-ray stack oi z ^ 2 BM/BX galaxies, 
Reddy & Steidel 2004), although this will be boosted by the lens. 
If it is an AGN, then we might expect the X-ray luminosity to be 
higher by a factor ~2 (Lehmer et al. 2005). In both cases, this has 
a negligible impact on our results. 

The most obvious feature of JEG3's environment is the ap- 
parent spatial offset between the radio galaxy itself and the bulk of 
galaxies at the group redshift of 0.649 (Fig. 6), coincident with the 
dominant source of X-ray emission. The galaxy's radio lobe mor- 
phology is suggestive of infall through a dense intracluster medium 
towards this higher density region. In Fig. 6 we represent the rel- 
ative velocities of galaxies about the central redshift, indicating 
objects with ±A«. Note that there are two other radio galaxies 



in this group; one ~0.5' to the west of the target radio galaxy, 
and the other (FR II source) several more arcminutes further to the 
west. Spectroscopic observations of these objects (Simpson et al. 
in prep) show that neither of them is associated with the structure 
surrounding JEG 3. These radio galaxies are not associated with the 
2 ~ 0.65 structure identified around the radio galaxy in this work 
(see Simpson et al. in prep). 

The radio galaxy appears to be associated with a small group 
interacting with the main body of the structure, thus it is not un- 
reasonable to associate the triggering of the radio galaxy with this 
interaction. This is a scenario which has been seen in other radio 
galaxy environments, at both high and low redshifts, and over a 
wide range of radio powers (Simpson & Rawlings 2002). For ex- 
ample, the high redshift radio galaxy TNJ1338-1942 {z = 4.1, 
Venemans et al. 2003) belongs to a 'proto-cluster', but is offset 
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from the centre of the overdensity. Similarly at low-redshift, the 
FRn source Cygnus A (z = 0.056) is involved in an cluster- 
cluster merger (e.g. Ledlow et al. 2005), with the radio-source off- 
set from the main overdensity (Owen et al. 1997). X-ray mapping 
of this cluster revealed hot (presumably) shocked gas between the 
two main peaks of the X-ray emission in the two sub-cluster units 
(Markevitch et al. 1999), consistent with the model of head-on clus- 
ter merging. These well studied radio galaxies are examples of the 
most luminous objects at their redshifts, whereas the luminosities 
of the radio galaxies in this work are modest in comparison. Trig- 
gering within interacting sub-cluster units seems to be responsible 
for a proportion of radio-galaxies, and this appears to have been 
happening at all epochs, and over a wide range of luminosities. 

The triggering of the radio galaxy could have occurred via 
galaxy-galaxy interactions in its immediate group, while the X-ray 
emission seems to be centred on the richer group in which gas is 
plausibly being heated by gravitational processes. However, there 
is a further complication to the Lx-Tx relationship, in that the out- 
come of the interaction between two sub-groups is a 'boosting' of 
the X-ray luminosity and temperature of the resulting cluster (Ran- 
dall, Sarazin & Ricker 2002). In general this will temporarily (i.e. 
for a few 100 Myr) enhance both Lx and Tx, before the cluster 
settles down to its equilibrium value. It is not clear what stage of 
interaction this environment is in. Given the relatively large off- 
set between the X-ray emission and the radio galaxy, we postu- 
late that this is in an early stage. In these environments it appears 
that we are witnessing radio galaxy activity and group-group merg- 
ing/interaction at the same time. Although this may hint at radio 
galaxy triggering in the cluster assembly phase, this does not imply 
that all clusters containing a radio galaxy are in an unrelaxed state 
- most relaxed clusters also contain a central radio source. This is 
interesting in that it appears that radio galaxy activity can be impor- 
tant in the thermodynamic history of the cluster environment over 
a wide range of scales (i.e. groups to massive clusters) and various 
stages of development. 

Following the assumption made in equation 4 of Miley (1980), 
and adopting an expansion speed of 0.03 — 0.15c (e.g. Scheuer 
1995), we estimate the jet power of JEG 3 to be of the order 0.5- 
1 X 10^® W, which is comparable to the X-ray luminosity of the 
cluster (the expected contribution to the X-ray luminosity from the 
radio galaxy itself is not significant). Although the similarity in 
power is likely coincidental, it is interesting, because the radio out- 
put of this low-power source appears to be sufficient to balance the 
radiative cooling of the gas. However, it is unlikely that the radio 
galaxy could keep up this balance over a long period: this mechani- 
cal injection will end long before the main merger between the two 
groups stops imparting energy to the system, given the expected 
timescales of radio activity (<10® yr, Parma et al. 2002) which will 
likely be the dominant event in the pre-equilibrium thermodynamic 
history of the cluster. Nevertheless, since the radio-galaxy appears 
to have been triggered early in the life-time of the cluster (the cen- 
tres of the radio and X-ray emission are separated by ~200kpc. 
Fig. 6), it could be important for providing significant energy input 
to the ICM, and therefore raising its entropy. 



5 SUMMARY 

We have presented multi-object (Low Dispersion Survey Spectro- 
graph 2) spectroscopy, combined with complimentary radio (VLA) 
and X-ray (XMM-Newton) observations of the Mpc-scale envi- 
ronments of a sample of low-power (Li,4ghz ^ 10^''WHz~^) 
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Figure 7. Spectra of components 'a' and 'b' of the radio galaxy in the field 
of JEG 3 - SXDF-iS- 170569. The inset grey scale images are 6" X 6" Sub- 
aiTj SuprimeCam thumbnails in BVRiz bands, clearly showing the dis- 
tinction between a red component 'a' and blue component 'b'. The spectra 
unambiguously identify 'b' with a galaxy at 2 = 1.847, strongly tensed by 
the foreground component 'a'. The emission lines at the very blue end of 
the spectrum 'b' are CiV, OlII] and Cm] in emission, securing the redshift. 
Fe III is seen in absorption at 2587A and 2600A. 



z ^ 0.5 galaxies in the Subaru-XMM Newton Deep Field (Simp- 
son et al. 2004). 

The radio galaxies targeted all appear to inhabit moderately 
rich groups, with remarkably similar X-ray properties, with JEG 3 
(VLA 0033 in Simpson et al. 2006) being the hottest and resid- 
ing in the richest environment. A statistical interpretation of the 
environment of JEG 2 (VLA 00 11 in Simpson et al. 2006) would 
suggest that it does not belong to a cluster, but this could be in 
part due to the spatial distribution of the galaxies - the other ev- 
idence that this is a relatively rich group is compelling: there is a 
relatively strong red sequence and extended X-ray emission seen at 
this location (suggesting a gravitationally bound structure), and the 
spectroscopic data (although dominated by small number statistics) 
suggests that this is indeed a group, or part of a much larger struc- 
ture, albeit one that is not relaxed. Also, three of the radio galaxies 
exhibit extended lobe morphology that suggest movement through 
a dense intracluster medium, again hinting that these systems are 
undergoing dynamical evolution. 

Our analysis concentrates on the environment of JEG 3, which 
presents us with a 'snap-shot' of radio-galaxy and cluster evolu- 
tion at 2 ~ 0.65. Overall, this is the richest environment we have 
studied, and the spectroscopic data reveal sub-structure: the radio 
galaxy belongs to a small group of galaxies which appears to be 
interacting (merging) with a larger group of galaxies. The X-ray 
emission for this system is associated with the richer, but radio- 
quiet group. We postulate that the radio emission has been trig- 
gered by a galaxy-galaxy merger within its local group. Although 
this may be providing some feedback to the ICM (indeed, the ra- 
dio output is capable of balancing the overall radiative cooling), 
the eventual group-group interaction will boost the overall X-ray 
luminosity and temperature. This process will take longer than the 
life-time of the radio-source, and so might be the dominating event 
in the thermodynamic fate of the environment. This scenario is not 
dissimilar to the situation with other radio galaxies at both high red- 
shift, and in the local Universe. Moreover, radio-triggering within 
sub-cluster units appears to be important for a proportion of the 
radio galaxy population over all epochs and luminosities. 
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